set term gif animate delay 1
set output "diffusion1d.gif"
set xrange [0:1]
set yrange [0:0.1]
set xlabel "x(m)"
set ylabel "T(K)"
do for[t=2:520:1]{
	plot "diffusion1dgif.txt" u 1:t w l lt 7 lw 3 t sprintf("Point source heat diffusion: t=%f s", (t-2)*0.00001)
}
set output
